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Abstract 

The viscous correction to thermal photon production at leading log order is calculated and in- 
tegrated over the space-time evolution of a hydrodynamic simulation of heavy-ion collisions. The 
resulting transverse momentum spectra and elliptic flow can be reliably calculated within a hydro- 
dynamic framework up to transverse momenta q± w 2.5 GeV and q± w 1.5 GeV respectively. A 
non- vanishing viscosity leads to a larger thermalization time when extracted from the experimen- 
tally measured inverse slope (T^q) of photon qj_ spectra. A precise, 0(20 MeV), measurement of 
photon Tgg- can place stringent bounds on tq and r]/s. 
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I. INTRODUCTION 



There is a general consensus that the matter produced at RHIC behaves as a near perfect 
fluid R. 0„e of the key flndiugs that led to this .nterpretatio,. is the stroug anisotrop.c flow 

of produced hadrons and its description by ideal hydrodynamic simulations [2| . Although it 
is too early to draw definitive conclusions it appears that the deviations from ideal hydro- 
dynamic behavior may be ascribed to dissipative effects. This has been suggested in many 
of the recent works on dissipative hydrodynamics jsj. 

Care must be taken when drawing conclusions about the viscosity of the early matter at 
RHIC based on hadronic observables alone. As hadrons interact strongly, they are sensitive 
to the later stages of the evolution leaving an ambiguity about whether the viscous effects 
seen in spectra stem from the hadronic or QGP phase. On the other hand, electromagnetic 
probes are emitted throughout the entire space-time evolution reaching the detector without 
any final state interactions and are not sensitive to the dynamics of freeze-out. 

As the next generation of experiments shift towards a precision study of the matter 
produced at RHIC it is imperative to have multiple observables constraining the medium 
properties inferred from the data. This work demonstrates that direct photons can be used 
to constrain the shear viscosity. It is also shown, that by neglecting the presence of viscosity, 
incorrect thermalization times will be extracted from experiments. 

There is a long history of photon calculations that we can't summarize here. Many of 
the works relevant to experiment jj] have relied on kinetic equilibrium and others have 
studied the effects of early momentum-space anisotropics. Only recently, however, has there 
been a measurement [g] that is precise enough to suggest the presence of an early hot stage 
of matter. 

II. PHOTON RATES WITH VISCOSITY 

In this section we show how the presence of viscosity modifies the photon spectra. In 



order to demonstrate the effect we look at the 2 — ?■ 2 processes in fig. [H 







'here are additional 



diagrams that contribute to the thermal photon rate at leading order [7|], which will not be 
examined in this leading log analysis. The emission rate of photons having momentum q 
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and energy Eg = qis 
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where the sum is over quark flavors and fa{Ea,Pa) is particle specie a's distribution function, 
which is not necessarily in equilibrium. We have used the upper (lower) sign for a final state 
boson (fermion). 
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FIG. 1: Feynman diagrams for the annihilation process (left) and the Compton process (right). 

For both the Compton and annihilation process the leading log behavior comes from when 
the exchanged momentum is soft (i.e. ~ gT). In this case the amplitude will be dominated 
by forward scattering. Following [Sj we can approximate the differential cross section as 

d(j 

Eg—^Eg(rtot{s)6\q~Pa), (2) 

where atot is the total cross section of either process. Performing the integration over pa we 
are left with 

W = jir^'^'^ ^ / ^'^^^^^^^^ ^ ^^^^^^^ ^ • 

By examining the above equation one can already see where the viscous correction will come 
into play. The expression for the photon rates is proportional to the distribution function 
of particle a (the quark). Therefore, if viscosity modifies the quark's distribution function 
it will modify the photon emission rates accordingly. Following {q] we take a one parameter 
ansatz for the viscous correction to the distribution function 

riPa) = fSiPa) + ^/o [1 ± /ol «5(M«.> , (4) 
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where /g is the particle's equihbrium distribution function and Ca is a constant determined 
in Appendix A and the notation (■ ■ ■ ) designates that the quantity in brackets should be 
symmetrized and made traceless. 

In general each distribution function in the above rate equation must be replaced by 
its viscous counterpart, fo + Sf, and then one must drop terms of 0{6p). A full analysis, 
which will be presented elsewhere, has found that the viscous corrections to the distribution 
functions occurring under the phase space integral lead to corrections to the coefficient under 
the log. Therefore, to leading log order, we neglect the viscous correction to and perform 
the above phase space integrals in the same manner as done in j^, [l^ where the diverging 
differential cross section is regulated using the re-summation technique^ of Braaten and 



Pisarski 



111. 



The final result of both the Compton and qq annihilation processes is 



" do 9 27r2 ■'"^^^ 



3.7388Eg 



where fa is the quark's off-equilibrium distribution function 



(5) 



riq) = fSiQ) + l-3|^/o [1 - /ol q'-Q'^d^.u^) . (6) 



III. SPECTRA IN ULTRA-RELATIVISTIC HEAVY-ION COLLISIONS 

In this section the above photon rates are integrated over the space-time evolution of 
the collision region determined by a 2+1 dimensional boost invariant viscous hydrodynamic 
model \l2\. The same bag model equation of state and Glauber model initial conditions of 

n 

[13l | are used since it was able to predict many hadronic observables. We have considered 
photon production from the QGP phase only. An impact parameter of 6 = 6.5 fm and 
Tj/s = 1/Att is used throughout. The one parameter in the rate equations, as, is evaluated 
at the scale fi = ttT from the two loop (3 function. Because we are using the leading-log 
rates the results below g± ~ 1 GeV are speculative. Above 1 GeV the expression under the 
log in eq. [5] remains larger than one. 



^ In principal one must also include the viscous correction when computing the resummed propagator. 
It turns out that these corrections can be taken into account by introducing a generalized momentum 
dependent thermal mass. This correction to the rates will not be enhanced by the logarithm and can 
therefore be neglected in this leading log approach. 
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There is one technicality regarding the viscous evolution model which must be discussed. 
We have chosen to use the ideal results for the evolution model for the viscous case as well. 
This amounts to neglecting the viscous corrections to the flow and temperature profiles 
which tend to be small, especially in the early stages of the evolution. Even though this 
approach is fundamentally inconsistent (energy-momentum is violated when converting from 
hydro to particles) the corrections are small. This procedure is convenient for two reasons. 
First, arbitrarily large gradients can be treated. And second, the final multiplicities will 
remain unchanged so we will not have to worry about modifying the initial conditions in 
order to re- fit hadronic observables at finite viscosity. The entire viscous effect, in this work, 
comes from the modification to the rates as outlined in the previous section. 



A previous work 



14| studied the emission of dileptons from a viscous medium taking 



into account the viscous correction to the underlying flow and temperature profiles. While 
the viscous correction to the hydrodynamic variables modified the yields the shape of the 
spectrum (as seen through Tgg) was largely undistorted. This is further motivation for 
neglecting the viscous corrections to flow. 

Fig. [2] shows the thermal photon transverse momentum spectra for the hydrodynamic 
model having starting times of tq = 0.2,0.6,1.0 fm/c. The ideal results are shown as 
solid lines. The dominant contribution at higher momentum comes from the early stages 
of the evolution where the temperature is highest. This leads to the expected increase in 
yields at higher momentum for earlier starting times. Shown as lines with symbols are the 
corresponding viscous results. A hardening of the transverse momentum spectra is seen. 



reminiscent of the known effect of viscosity on single particle spectra [15j. These results 
should not come as a surprise. At earlier times the collision geometry has larger gradients 
and the underlying quasi-particles are furthest from local thermal equilibrium. Viscosity 
introduces a power law correction to the spectra (reminiscent of particle production in 
perturbative QCD). As the system evolves, thermalization occurs by transferring momentum 
by bremsstrahlung or collisions to softer modes until the spectra eventually become thermal. 
Since the relaxation time grows with energy incomplete thermalization enhances the quark 
distribution at high qj-. The harder distribution of quarks leads to a harder spectrum of 
photons. 

The dominate contribution to the transverse momentum spectra at high q± is from the 
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FIG. 2: q± spectra of thermal photons from the QGP (top) and the ratio of the viscous correction 
over the ideal result (bottom). The band in the lower figure indicates where the q± spectra can no 
longer be reliably calculated. 

earliest times when the medium is hottest and the gradients are largest^. This leads to a 
larger viscous correction, and the subsequent breakdown of the hydrodynamic description, 
at large q±. The lower plot of fig. [2] shows the ratio of the viscous correction to the ideal 
result. The correction becomes of order one at q± ^2.5 GeV. 
Fig. |3] shows the photon elliptic flow defined as 

/ c/0cos(20) dN + 6dN J rf0cos(20) dN + 5dN J d<p 5dN J rf0cos(20) dN 



V2 = 



(7) 



Jd(pdN + 5dN J d<p dN (/ d<p dNf 

where dN = dN/ d^q is the space-time integrated ideal photon spectra and 6dN = 5dN/ d^q 
is the viscous correction. In the rightmost expression the denominator has been expanded 



The viscous correction is proportional to 9^^'' ~ l/r at early times 
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FIG. 3: Elliptic flow of thermal photons from the QGP. 



in order to keep terms up to 0{6f) only. The solid red line shows the ideal result. As 
explained in 16|, the total photon V2 follows from a weighted average of the flow over the 
proper time of the evolution. A linear rise in V2, as expected from hydro, is observed. The 
suppression at large q± is due to the early non-flowing phase, which dominates the yields at 
large momentum. Including viscosity, shown as the line with symbols, has a large effect on 
the f 2. The ideal result is suppressed by a factor as large as ~ 2. In contrast to the hadronic 
V2, where the largest effect of viscosity is at high we find large corrections at all qj_ in 
the case of photons. 

The solid magenta line in Fig. [3] shows the viscous result using the expansion in the 
rightmost expression of eq. [71 When the two results disagree this observable can no longer 
be reliably calculated. This happens when q± ^ 1.5 GeV as shown by the solid band. 



IV. DISCUSSION 



The sensitivity of the spectra to the shear viscosity makes the extraction of the thermal- 
ization time more subtle. It also makes it possible to extract rj/s from the data. In order 
to demonstrate the process, the inverse slope of the photon spectra is computed by a fit to 
l/q_idN/dq± cx exp(— g_|_/Tgg) in the momentum region'^ 1.5 < q± (GeV) < 2.5. In fig. H] the 

^ Realistic model calculations find that the QGP contribution dominates the yields in this kinematic 
region. 
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FIG. 4: Effective temperature of final photon spectra for starting times of tq = 0.2,0.6,1.0 and 
1.4. Lines are shown to guide the eye. The sohd band is the measured slope from minimum bias 
data ^. 

effective temperature T^g is plotted versus the thermalization time tq for rj/s = 0,1/Att and 
2/ Ait. Both viscosity and earlier thermalization cause the effective temperature to increase 
in a non-trivial way. 

We should clarify to the reader that what we call thermalization time is really the hy- 
drodynamic starting time. In this work any production before tq has been neglected. If one 
included non-equilibrium production from the early evolution prior to tq we would expect 
to see the strong dependence on tq to be reduced. 

From fig. HI it is clear that fitting the data with ideal hydrodynamic simulations will result 
in the extraction of earlier thermalization times. For example, the same T^g is observed for 
an ideal (viscous) evolution starting at tq =0.6 (1.0) fm/c. 

A precise measurement of the inverse slope could constrain a combination of tq and t]/s. 
In order to discern between the ideal and viscous results shown here a measurement must pin 
down Tgg to within 20 MeV. The band in fig. H] shows the experimentally measured slope |6|, 
including both systematic and statistical errors, in order to demonstrate the current quality 
of the data. 

There are a number of caveats which must be discussed before a fair comparison can be 
made with data. The largest uncertainty comes from using the leading log results. Exam- 
ining the full leading order results from [3] we estimate that the leading log contribution 
used here comprises about fifty percent of the thermal QGP photon yields in the kinematic 
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regions of interest. Therefore, even if the viscous correction to the remaining leading or- 
der (non leading log) results turned out to be neghgible, one would still expect the above 
conclusions to hold at a qualitative level. 

V. CONCLUSIONS 

The viscous correction to thermal photons at leading log order is computed. The rates 
are then integrated over the space-time history of a hydrodynamic simulation. The viscosity 
sets a bound on where a hydrodynamic description is reliable {q± ~ 2.5 GeV for spectra 
and q± ~ 1.5 GeV for elliptic flow). We have shown that viscosity increases T^q and by 
neglecting its presence wrong conclusions about the thermalization time will be reached. 
Our model calculation has shown that the photon spectra can place stringent constraints on 
both To and r]/s. 
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Appendix A: Strength of the off-equilibrium quark distribution 

In this appendix we make an estimate for Cg. From the definition of the viscous part of 
the stress-energy tensor 




(Al) 



one can obtain a relationship between the viscosity of the medium and the constants^ Cg 



and C, 




(A2) 



* The sum is now restricted over the two species a = q (quark), (gluon) of the QGP. 
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In the above expression s = Sq + Sg is the total entropy and Ng = ^^f^^i^^ ~ 1.031 and 
Ng = ^^^3^ ~ 0.958 are two constants that come from the phase space integrals over the 
quark and gluon distribution functions. The above expression is simply the statement that 
Vtot = Vq ~^ Vg- order to solve for Cq we need to estimate the relative viscosity from 
the quark and gluon phases. Following {q] we take rjg ^ 1. 70 Nf rjg. This factor can be 
understood from the relative strength of the three cross sections relevant to viscosity at 
leading log order; gg — )• gg:qg — )■ qg:qq — > gg = | : 2 : |. 

For Nf = 2 we find Cq ~ l.Srj/s. The larger coefficient for quarks intuitively makes sense. 
The quarks must be even further out of equilibrium in order to compensate for the quicker 
relaxation of the gluons. 
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